Interactive MAP:

Click the point to popup information

leaflet() %>% 
  #  addMiniMap(position = "bottomleft") %>%
  addProviderTiles("Esri.WorldImagery") %>% 
  addPolygons(data = area1, color = "white", weight = 1, fillOpacity = 0) %>% 
  addCircles(data = loc_trees, color = "yellow", opacity = 0.9,
             popup = paste("Specie:", loc_trees@data$species, "<br>",
                           "Kind of damage:", loc_trees@data$tipo_de_dano, "<br>",
                           "Succesional stage:", loc_trees@data$succ, "<br>",
                           "Density:", loc_trees@data$densidade))

Position and Density:

ggplot(df, aes(x=x, y=y))+
  geom_point(aes(col = parcela_bl), size = 3, alpha = 0.5)+
  scale_color_manual(values = eqm)+
  ggtitle("Fall Trees Location")+
  theme_light()+
  coord_fixed()
## Warning: Removed 5 rows containing missing values (geom_point).

ggplot(df, aes(x=x, y=y))+
  stat_density2d(geom="tile", aes(fill = ..density..), contour = FALSE) + 
  stat_density_2d(col = "black")+
  scale_fill_viridis_c()+
  ggtitle("Fall Trees Density")+
  theme_light()+
  coord_fixed()
## Warning: Removed 5 rows containing non-finite values (stat_density2d).
## Warning: Removed 5 rows containing non-finite values (stat_density2d).

Climax and Pionner trees:

int = ggplot(df, aes(x=x, y=y))+
  geom_point(aes(col = succ), size = 3, alpha = 0.5)+
  scale_color_manual(values = eqm2)+
  ggtitle("Fall Trees Succecional Condition")+
  theme_light()+
  coord_fixed()

ggplotly(int)